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FORTRAN PROGRAM FOR CALCULATING AERODYNAMIC FORCES FROM 
PRESSURE OR VELOCITY DISTRIBUTIONS ON BLADE SECTIONS 

by William D. McNally 
Lewis Research Center 

SUMMARY 

A FORTRAN IV program is presented which calculates aerodynamic forces on tur- 
bomachinery blade sections from distributions of pressure or velocity along the sur- 
faces. Forces due to pressure on either single blade sections or tandem blade sections 
may be calculated. 

The program is designed to be used in conjunction with other existing programs 
which calculate blade surface and internal velocity distributions, streamlines, and 
boundary layers. Together these programs permit the analytical study of the flow 
through turbomachinery blading. 

Input to this program includes blade surface coordinates, surface distributions of 
pressure or velocity, and several overall flow parameters. Spline curves are fitted to 
the surface coordinates to obtain surface angles; or these may be given as input. Pres- 
sure or velocity distributions are integrated around the blades using the trapezoidal rule 
to obtain components of force in the input coordinate system directions. Meridional and 
tangential forces, and lift and drag forces (due only to pressure) are computed from the 
integrated force components. For tandem blade sections, ratios of rear blade segment 
to front blade segment forces are also given. 

Care must be exercised in the placement of points near blade leading and trailing 
edges. The points in these areas have high values of surface angle with corresponding 
large tangents used in the integrations. For this reason, input points must be more 
closely spaced in these regions. 

This report includes a listing of the FORTRAN IV computer program, with an ex- 
planation of the input required and the output generated. Three tandem blades with dif- 
ferent ratios of camber between the blade segments are analyzed with the program. A 
discussion is included concerning the placement of points on blade surfaces in order to 
obtain accurate results with the program. Running times are about 1/10 minute on IBM 
7094 equipment. All parts of the program are in general FORTRAN IV code, and could 
be easily transferred to other IBM equipment. 



INTRODUCTION 


In the design, analysis, and comparison of blades for turbomachinery, it is desir- 
able for a number of reasons to have knowledge of the aerodynamic forces acting on 
blade sections: 

(1) Lift forces, or tangential forces, are directly related to the energy addition or 
work input. 

(2) Blade surface forces due to pressure give a measure of each blade's "loading, " 
where loading is thought of in terms of area under blade surface pressure profiles or 
total force on the blade section. 

(3) For tandem blade sections, tangential forces provide a measure of "work-split" 
between the front and rear segments of the blades. 

(4) Blade forces are used for the calculation of aerodynamic stresses and moments 
in the mechanical design process. 

Forces on turbomachinery blading are of two kinds: pressure forces (mainly lift) 
and skin friction forces (mainly drag). In ideal flow, skin friction is neglected and the 
drag force is zero. So the only force is lift due to pressure. In real flow, skin friction 
exists and the drag is not zero. The component of drag due to skin friction can be com- 
puted with some degree of accuracy from boundary layer theory. Besides the skin fric- 
tion drag, however, there is another type of drag called form drag or induced drag or 
pressure drag. This drag component due to pressure arises from the boundary layer 
displacing the free stream from the blade and modifying the ideal flow pressure distri- 
bution. Pressure drag is difficult to calculate (see refs. 1 and 2), especially from inte- 
gration of real flow pressure profiles. On the other hand, pressure forces in the lift 
direction can be obtained from integration of ideal flow pressure profiles. 

There are several methods available for obtaining forces on blade sections. One 
method is through the use of vector diagrams across the blade section with the assist- 
ance of empirical relations for fluid turning and loss. An alternate method, used more 
in the analysis of today's highly loaded blading which has not yet been tested, is to ob- 
tain forces analytically from the integration of blade surface pressure profiles. 

This report describes a computer program developed to calculate aerodynamic 
forces on turbomachinery blading from given surface distributions of pressure or ve- 
locity. It is recommended primarily for use with ideal flow pressure distributions 
where forces in the drag direction can be neglected. 

The force program described herein is particularly useful in conjunction with exist- 
ing computer programs discussed in references 3 to 6. These programs calculate sur- 
face pressure distributions, free-stream and internal velocity distributions, stream- 
lines, and boundary layers. Combining the results of all these programs allows the 
analytical study of the performance and two-dimensional flow through turbomachinery 
blading before it is studied experimentally in the laboratory. 
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This report presents a brief derivation of equations, a listing of the program, and a 
description of its input and output. Input to the program is easily obtained from the out- 
put of other programs for computing pressure or velocity distributions. Geometrical in- 
put consists of X-Y coordinates of blade surfaces in any two mutually perpendicular di- 
rections. Surface angles may be given as input or may be computed in the program from 
spline curves. Output consists of spline curve fits of surface coordinates and integra- 
tions of pressure along each blade surface, as well as integrated blade forces in desired 
directions. The report also includes numerical examples to illustrate typical input 
values and the form in which output is given. Internal variables used in the program are 
defined only if they pertain to input. 


GENERAL DESCRIPTION OF PROGRAM 


The program computes blade forces due to pressure distribution for either a single 
blade section or a tandem blade section (with two blade segments). Each blade section 



M 



Figure 1. - Surfaces on which blade surface velocities 
are known. 


Y 




(b) X axis through leading edge or trailing edge circle centers. 



(c) X axis tangent to pressure surface leading edge and trailing edge radii. 
Figure 2. - Typical variations in placement of X-Y input axes. 
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(i.e. , each "cut" through a single or tandem blade) requires a separate set of input. 

Input may be given on any plane on which blade pressures or velocities are known. 
Usually these are obtained on a cylindrical surface (Z -6 coordinates) or on an axisym- 
metric surface of revolution (M-0 coordinates) (see fig. 1). These unwrapped surfaces 
are the input planes for this program. 

Blade geometry is defined by X-Y coordinates of each blade surface. The X and 

Y axes can be in any two mutually perpendicular directions, and are related to the Z -9 

or M-0 directions by the input angle (see fig. 2). As figure 2 shows, the X and 

Y axes can be placed on the blade in any of the standard positions. 

The second portion of the input, a description of the flow, is given by means of 
either blade surface static pressure or blade surface velocity distributions. Figure 3 


R 



indicates the blade surface velocity components in the input planes: V z and V 0 on the 
cylindrical blade-to-blade plane, and V M and Vq on the general surface of revolution 
blade-to-blade plane. 

Input is given on one blade surface at a time in the order indicated by figure 4. For 
single segment blades, only surfaces 1 and 2 are used. For a tandem blade, four sur- 
faces are used and the front blade segment is entered before the rear blade segment. 

The integration scheme requires blade surface angles (with respect to the X axis) 
at the X-Y input points (see fig. 5). These may either be given by the user as input, or 
be computed by the program by means of spline curve fits of the input geometry (refs. 7 
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Figured -Numbering of Wade surfaces for input. 



Figure 5. - Input parameters on typical surface. 


and 8). The user may also give blade angles at a portion of the input points. In this 
case the program uses these input values where they are given, and uses computed 
angles from the spline curves at all the other points. In all cases spline fits of the blade 
surfaces are made; and slopes, curvatures, and angles are printed in the output. 

Both blade surface pressures and velocities are printed as output. The relation be- 
tween the two that is used in the program is Bernoulli's equation 


y /pV p \ ( y~ 1}/ y , y 2 _ y /p’\ 
y-lVp’A p’/ 2 y - 1 \p' / 


(1) 


Force components are obtained initially for each blade surface by integrating pres- 
sure distributions in the X and Y directions. On the suction surface the force in the 
X direction is given by 
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X, s 


Z 


P„ AY, 


Z P S 


cn AX_ 
s s 


( 2 ) 


Using the average surface angles over dX increments, equation (2) becomes at each 
point i 




(3) 


The force in the Y direction on the suction surface is 


F 


Y,s 


Z'- P s> ix s 


(4) 


which becomes 




-( Ps ’ il 2 Ps,1 - 1 ) (X s,i- X s,i-l ) 


On the pressure surface the force in the X direction is given by 


(5) 


F x,p = Z ( - p P )AY P := Z ( - I> P ),an “p ax p 


( 6 ) 


Again, using average surface angles, equation (6) becomes at each point 


F X, p, j 
j=l 


i-1 

F x, P ,j 

i=i 



The force in the Y direction on the pressure surface is 


(7) 


r Y,p=Z P p 4X p 


( 8 ) 
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which becomes 




Y,p, j 


+ 





p,i-l ) 


(9) 


After force components are obtained in the X and Y directions on each blade surface, 
they are rotated to the M and B directions (see fig. 2) as follows 


F M = F X cos a MX ~ F Y sin “MX 
F 0 = F X sin “MX + F Y cos “MX 


( 10 ) 


Finally force components are rotated to the lift-drag (L-D) directions (perpendicular and 

parallel to the direction of mean velocity V at the mean flow angle B , fig. 6) 

° 1 mean 7 ° 

as follows 


F D = F M cos ^mean + F e sin ^mean " 
F L = - F M sin ^mean + F 0 cos ^mea n> 


(ID 



M or Z 


R0 




Inlet 


Mean velocity and flow angle across 
blade section 


Outlet 



Figure 6. - Mean velocity vector and flow direction notation. 
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The mean flow angle across the blade section P mean 
fig. 6) 


is computed by the formula (see 


P 


mean 



+ ^out, 8 \ 
+ V out, Z / 


( 12 ) 


After the calculation of force components (in various directions) on the blade sur- 
faces, the components are summed for pressure and suction surfaces to obtain total 
forces in certain directions on each of the blade segments . In the X- Y directions these 
forces are obtained as follows 


(13) 

Similar expressions are used in the M-0 and L-D directions. 

After X-Y, M-0, and L-D forces are computed for each blade segment, total 
forces and angles with respect to the M axis are computed for the segments . At this 
point the program is finished for a single segment blade section. 

For tandem blades the program continues and adds together forces for both blade 
segments in all the principal directions. Ratios are also given for forces on the rear 
blade segment to forces in similar directions on the front blade segment. 

For an ideal flow pressure distribution, the final drag force F D for either a single 
segment blade or a tandem segment blade should be very close to zero. Inaccuracies in 
the input ideal flow pressure distribution and numerical inaccuracies in the program pre- 
vent Fj-j from being exactly zero. For a tandem blade , since the L and D directions 
are perpendicular and parallel to the direction of mean velocity across the overall blade 
section, F^ for either of the individual segments will not usually be zero. But F^ 
should be near zero when the segment contributions are summed for the overall blade 
section. 

Output from the program consists of the following: (1) printing of all the input, 

(2) spline curve fit data for each surface, (3) summations from point to point of the sur- 
face pressure distributions to obtain force components, and (4) summed forces in all the 
principal directions. 

The program is written entirely in FORTRAN IV and is run at NASA Lewis on the 
IBM 7094-7044 direct coupled system with a 32 768 word core (77777^). The total pro- 
gram storage requirement is 22436^ of which 6335^ is used in the storage of vari- 
ables. The program runs in 1/10 minute on IBM 7094 equipment. 

At Lewis, the program is currently being used in conjunction with the programs of 


F X ~ F X, s + F X, p 


F Y - F Y,s + F Y,p 
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references 3 to 5 to compute forces and evaluate work input and work splits on compara- 
tive blade designs. Ideal flow velocity distributions are being used as input. 


NUMERICAL EXAMPLES 

Three blade sections are included as numerical examples to illustrate the use of the 
program. All three are tandem blade sections with varying ratios of camber of the rear 
segment to camber of the front segment. These blades were generated by means of the 
program of reference 9. Table I contains some of the characteristic parameters of these 
blades. As table I indicates, the camber ratios vary from 1:1 to 3:1 (i.e., camber of 
rear blade segment three times camber of front blade segment). For all three blades 
the chord length of the rear blade segment is equal to the chord length of the front blade 


TABLE I. - CHARACTERISTIC PARAMETERS OF THREE TANDEM 
BLADES FOR NUMERICAL EXAMPLE 


Tandem 

blade 

number 

Camber 

ratio 

Front blade 
camber, 
deg 

Rear blade 
camber, 
deg 

Total chord 

Total 

camber, 

deg 

Inlet 

blade angle, 
deg 

Solidity, 

chord/spacing 

ft 

m 


in 

39.2 

39.2 

0.3333 


58.8 

46.25 

■m 


SHI 

26.2 

52.4 

.3333 


58.8 

46.25 



K8fl 

19.6 

58.8 

.3333 


58.8 

46.25 

m 


segment, and the overlap is 20 percent of the chord of the first blade segment. All three 
blade sections have a converging channel between blade segments, with the openings at 
the front and rear of the channel equal to 14 and 10 percent, respectively, of the front 
blade chord. 

Figure 7 shows the three blade geometries superimposed. Figure 8 contains the 
three surface velocity distributions obtained by means of the program of reference 4. 

The input required for obtaining forces on blade number 2 is listed for illustration in 
table II; the variables are explained in the section Input Variables. Table III lists the 
important blade forces on each of the three tandem blades obtained by means of the pro- 
gram described herein. Lift, drag, and total force are also given for the entire tandem 
blade sections. Ideally, the drags would be zero, but these small components of drag do 
not noticeably affect the total forces on the blades. One of the principal points to notice 
on the table is the shift of forces from front to rear blade segment. The magnitude of 
these forces provides a measure of the work split F 0 and loading split F tot between 
the two blade segments. 
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Surface velocity, ft/sec Distance in R0 direction, ft 
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TABLE H. - INPUT FOR TANDEM BLADE NUMBER 2 
OF NUMERICAL EXAMPLES 


T AN OEM OCA BLADE C2/C1=1.0 PH2/PHl=2.0 G/Ci=.l L/C1=.2 F=1.4 


GAM 

R 

PTZ 

TTZ 


1.400 

1 716.48 

2080.00 

520.00 


VR 1 

VRO 

BETAI 

BET AO 

ALPHMX 

351. 291 00 

249.63800 

43.250 

-1.000 

0. 

KP V 

KSIRF 

NSURF1 

NSURF2 

NSURF3 

2 

4 

27 

36 

28 


SUCTION SURFACE - FRONT BLADE 


XI 

Y 1 

VI 

ANG1 

0. 

0. 

389.0000 

-0. 

0.000400 

0.001150 

389.7Q00 

-0. 

0.001000 

0.001920 

39G.5Q00 

-0. 

0.002000 

0.003240 

392.3000 

-0. 

0.005000 

0.007050 

397.1000 

-0. 

0.0 15164 

0.019214 

411.2QOO 

-0. 

0.025274 

0.029916 

422.7Q00 

-0. 

0.035384 

0.039514 

430.0000 

-0. 

0.045493 

0.048186 

433.4000 

-0. 

0.055603 

0.056028 

432.0000 

-0. 

0.065713 

0.063142 

429.0000 

-0. 

0.075822 

0.069577 

424.0000 

-0. 

0.085932 

0.075391 

417.7000 

-0. 

0.096041 

0.080644 

410.400Q 

-0. 

0.106150 

0.085368 

401.0000 

-0. 

0.1 16260 

0.089584 

390.5000 

-0. 

0.126370 

0.093328 

377.9Q00 

-0. 

0.136480 

0.096622 

365. 2Q 00 

-0. 

0.143130 

0.098534 

353.7000 

-0. 

0.145770 

C. 100268 

341 .2000 

-0. 

0.152000 

0.100810 

336.9000 

-0. 

0.153100 

0.101078 

335.0000 

-0. 

0.1 54000 

0.101270 

333.0000 

-0. 

0. 1 55000 

0.101500 

331.0Q00 

-0. 

0.155800 

0.101300 

329.3Q00 

-0. 

0.1 56200 

C. 100930 

328.5000 

-0. 

0.156420 

0.100088 

328.0000 

-0. 

PRE SSURE 

SURFACE - FRONT 

BLADE 


X2 

Y2 

V2 

ANG2 

0 . 

0. 

286.0000 

-0. 

0.000400 

-0.001090 

285.8000 

-0. 

0.0C1000 

-0.001600 

285. 5Q 00 

-0. 

0.002000 

-0.001800 

285.0000 

-0. 

0.003000 

-0.001470 

284.6000 

-0. 

O.0C4000 

-0.000670 

284.2000 

-0. 

0.005000 

0.000150 

283.9000 

-0. 

0.006000 

0.000940 

283.4000 

-0. 


NSURF4 

36 
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TABLE H. - Continued. INPUT FOR TANDEM BLADE NUMBER 2 


OF NUMERICAL EXAMPLES 


O.CC8O0O 

C. 002580 

282.9000 

-0. 

0.0 15164 

0.008266 

281.0000 

-0. 

0.025274 

0.016086 

278. 1000 

-0. 

0.035384 

0.023694 

273.4000 

-0. 

0.045493 

C. 031025 

268.2000 

-0. 

0.055603 

0.038250 

263.4000 

-0. 

0.065713 

0.045199 

258.5000 

-0. 

0.0 75822 

0.051956 

253.0000 

-0. 

0.080877 

0.055250 

249.6000 

-0. 

0-085932 

0.058522 

246.6000 

-0. 

0.090987 

0.061731 

243.3Q00 

-0. 

0.096041 

0.064898 

239.8000 

-0. 

0.101100 

0.068021 

236.9000 

-0. 

0.106150 

0.071102 

234.7Q00 

-0. 

0.111210 

0.074120 

232.4000 

-0. 

0.114000 

0.075750 

232.0000 

-0. 

0.1 16260 

0.077095 

232.2000 

-0. 

0. 121320 

C. 080028 

234.3000 

-0. 

0.126370 

0.082918 

239.9000 

-0. 

0.131430 

0.085744 

247.4000 

-0. 

0.136480 

0.088549 

258.5000 

-0. 

0.143130 

0.092210 

278.0000 

-0. 

0.149770 

0.095752 

300.0Q00 

-0. 

0.153000 

0.097480 

312.9Q00 

-0. 

0.155000 

0.098440 

321.1000 

-0. 

0.155 800 

0.098850 

324.8000 

-0. 

0.156200 

0.099240 

326.8000 

-0. 

0.1 56 420 

0.100088 

328.0000 

-0. 

SUCTION SURFACE - REAR 

BLADE 


X3 

Y3 

V3 

ANG3 

0.136480 

0.057651 

218.0000 

-0. 

0.136800 

0.058680 

219.8000 

-0. 

0. 137400 

C. 059400 

223.0000 

-0. 

0.138000 

0.060000 

226.2Q00 

-0. 

0.1 39 800 

0.061838 

235.0000 

-0. 

0.146450 

0.068170 

263.1100 

-0. 

0.153100 

0.073822 

287.4000 

-0. 

0.161 460 

0.080198 

312.0Q00 

-0. 

0.1 71530 

0.086864 

333.2000 

-0. 

0.181600 

0.092556 

352.0000 

-0. 

0.151670 

0.097393 

369.1000 

-0. 

0.201740 

0.101454 

378.4Q00 

-0. 

0.2 11810 

0.104782 

384.0000 

-0. 

0.221880 

0.107427 

386.0000 

-0. 

0.2 31 550 

0.109435 

385.5000 

-0. 

0.242030 

C. 110821 

381.0000 

-0. 

0.252100 

0.111592 

373.3000 

-0. 

0.262170 

0.111754 

363.5000 

-0. 

0.272240 

0.111320 

353.0000 

-0. 

0.282310 

0.1 10292 

341.6000 

-0. 

0.292380 

0.108662 

324.2000 

-0. 

C.3C2450 

0.106382 

300.9000 

-0. 
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TABLE H. - Concluded. INPUT FOR TANDEM BLADE NUMBER 2 


OF NUMERICAL EXAMPLES 


0.312520 

0.103436 

269.8000 

-0. 

0.3 15000 

0.102600 

259.5000 

-0. 

0.3 16000 

0.102220 

254.8000 

-0. 

0.3 16600 

0-102000 

252.0000 

-0. 

0.317200 

0.101600 

249.0000 

-0. 

0.317560 

0.100725 

247.4000 

-0. 

PRESSURE 

SURFACE - REAR 

BLADE 


X4 

Y4 

V4 

ANG4 

0.136580 

0.057651 

480.0000 

-0. 

0.136800 

0.056680 

473.0000 

-0. 

0. 1 37400 

0.056120 

458.0000 

-0. 

0.138000 

0.055870 

442.5000 

-0. 

0.139000 

0.055900 

418.4000 

-0. 

0.140000 

0.056430 

398.2000 

-0. 

0.141000 

0.057110 

379.2000 

-0. 

0.143000 

0.058400 

343.5000 

-0. 

0.146400 

0.060540 

290.0000 

-0. 

0.148000 

0.061490 

2 71.4000 

-0. 

0.150000 

0.062670 

257.0000 

-0. 

0.151600 

0.063610 

249.5000 

-0. 

0.153100 

0.064450 

245.6000 

-0. 

0.155000 

0.065560 

241.7000 

-0. 

0.161460 

0.068999 

233.4000 

-0. 

0.171530 

0.074014 

222.1000 

-0. 

0.181600 

0.078540 

215.5000 

-0. 

0.191670 

0.082578 

211.0000 

-0. 

0.201740 

0.086169 

2 07.0000 

-0. 

0.211810 

0.089335 

2 03. 7 Q 00 

-0. 

0.221880 

C. 092098 

200.8Q00 

-0. 

0.231950 

0.094456 

199.8000 

-0. 

0.242030 

C. 09641 1 

199.5000 

-0. 

0.252100 

0.097984 

199. 7000 

-0. 

0.262170 

0.099174 

200. 8Q00 

-0. 

0.272240 

C. 100002 

202. 7000 

-0. 

0.282310 

0.100491 

205.6000 

-0. 

0.292380 

C. 100576 

210.8000 

-0. 

0.302450 

C.l 00300 

219.9000 

-0. 

0.308000 

0.100080 

230.0000 

-0. 

0.3 12520 

0.099662 

242.3000 

-0. 

0.3 15000 

0.099400 

251.0000 

-0. 

0.316000 

0.099290 

255.0000 

-0. 

0.316600 

0.099370 

257. 7000 

-0. 

0.317200 

C. 099790 

260. 1000 

-0. 

0.317560 

0.100725 

261.3000 

-0. 
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Figure 9. - Input form. 


INPUT 

Figure 9 shows the input variables as they are punched on data cards. The first 
input card is for a title which identifies the data deck and is printed on the output. The 
user may type whatever information he wishes in any of the columns of this card. The 
remaining cards are for input data; the input variables are defined in the next section. 


Input Variables 

GAM specific heat ratio, y 

R gas constant, R, (ft)(lbf)/(slug)(°R); j/(kg)(K) 

PTZ inlet or upstream relative total pressure, P', (lbf)/(ft 2 ); 

N/(m) 2 
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TTZ 

VRI 

VRO 

BETAI 

BETAO 


ALPHMX 


KPV 

KSURF 


NSURF1(NSURF2, 
NSURF3 , NSURF 4) 

X1(X2,X3,X4) 
Y1(Y2, Y3, Y4) 


P1(P2, P3, P4) 
VI (V2, V3, V4) 


inlet or upstream relative total temperature, T*, °R; K 

inlet or upstream flow velocity relative to blade row, Vj , 
ft/sec; m/sec (figs. 4 and 6) 

outlet or downstream flow velocity relative to blade row, 
V out’ ft / sec > m / sec (figs. 4 and 6) 

inlet or upstream relative flow angle, /3^ n , deg (figs. 4 and 6) 

outlet or downstream relative flow angle, /3 ., deg (figs. 4 

and 6) BETAO is positive if rotation from the M or 
Z axis to V out is counterclockwise; negative if rotation 
is clockwise. 

angle between meridional axis and X axis of input, 
deg (fig. 2) -90° < < ' 90° > ALPHMX is positive if 

rotation from the M or Z axis to the X axis is counter- 
clockwise; negative if rotation is clockwise. 

integer (1 or 2) indicating whether pressures (1) or velocities 
(2) are used to input flow conditions on the blade surfaces 

integer (2 or 4) indicating whether input is given for a single 
segment blade section (two surfaces) or for a tandem blade 
section (four surfaces) 

integer number of input points on surface 1 (surfaces 2 to 4) 

If KSURF = 2, NSURF3 and NSURF 4 should be zero. 

array of X coordinates of input points on surface 1 (sur- 
faces 2 to 4), X, ft; m (figs. 2 and 5) 

array of Y coordinates of input points on surface 1 (sur- 
faces 2 to 4), Y, ft; m (figs. 2 and 5) The X and Y co- 
ordinates on all surfaces are given with respect to the same 
X-Y origin. There is not a separate origin for each blade 
segment (fig. 4). 

array of static pressures at the input points on surface 1 
(surfaces 2 to 4), P, lbf/ft^; N/m^ (fig. 5) 

array of relative surface or free-stream velocities at the 
input points on surface 1 (surfaces 2 to 4), V, ft/sec; 
m/sec (see fig. 5) Either pressure P or velocity V is 
used to input flow conditions on the blade surfaces. KPV 
indicates which of the two is used in a set of input data. 
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ANG1(ANG2, ANG3, ANG4) array of surface angles with respect to the X axis at the 

input points on surface 1 (surfaces 2 to 4), a, deg (see 
fig. 5) The angle a is positive if rotation from the X 
axis to the tangent to the surface is counterclockwise; 
negative if rotation is clockwise. Angles may be given at 
all points, at no points, or at only some of the points. If 
any angles are given, spline curve results will be overriden 
at the points where angles are given as input. 

All the information for a single input point (X, Y, P or V, and a) is given on a single 
data card (see fig. 9). For a single blade (KSURF = 2), NSURF3 and NSURF4 are set to 
zero, and no cards are given for blade surfaces 3 and 4 (X3, Y3, etc. , and X4, Y4, etc.) 

Placement of Points for Accurate Results 

In order to obtain accurate values for force, particularly in the drag direction, sev- 
eral factors have to be kept in mind. First, consider the trapezoidal integration scheme. 
It approximates areas under a curve of the function to be integrated by areas of trape- 
zoids. In order to achieve accuracy, sufficient points are needed in regions where the 
slope of the function is changing rapidly. The function to be integrated in this program 
depends upon pressure and surface angle, and these change most rapidly around the lead- 
ing and trailing edges of the blades. The computer program which generates the ideal 
flow pressure or velocity distribution may not give enough points in these regions. If 
that is the case, plots must be made of the blade leading and trailing edges and the pres- 
sure or velocity distributions in these regions in order to interpolate data for input. 

This procedure had to be followed for the blades used in the numerical examples of this 
report. Figure 10 shows this plot for the front blade segment of blade number 2 of the 
numerical examples. The circled points indicate where surface velocities were available 
from the ideal flow program of reference 4. Points indicated by triangles were added for 
use in this program. 

A second factor is that the velocities plotted in figure 10 are in error very close to 
the leading edge of the blade. The velocity distributions do not include the deceleration 
to zero at the stagnation point, nor the peak velocities occurring due to rapid accelera- 
tions around the leading edge radius. Figure 10 merely contains enlarged portions of the 
plot of figure 8 which represents the output of the programs of reference 4. The very 
local peaks about the leading and trailing edges of this blade do not affect the overall 
loading significantly, and can be neglected. If need be, they can be calculated using ref- 
erence 5 and included in this analysis. 

A third point is that the force in the X direction (generally close to the drag force 


17 





Distance in M direction, ft 


Figure 10. -Leading and trailing edge surfaces and velocities for front blade segment of tan- 
dem blade number 2. 


direction) is very sensitive to the blade surface angles at the input points, particularly 
around the leading and trailing edges. For example, changing the number of points 
around a leading edge can alter the drag force as much as 25 percent while the lift force 
remains essentially unchanged. This is why the leading edge must be described in detail 
as it has been on figure 10. In general, 25 to 35 points has been found to be adequate for 
describing most blade surfaces, with about one-third of these concentrated near the lead- 
ing and trailing edges. 

Finally, attention must be given to the spline curve fits of the blade geometry. If 
blade surface angles are not supplied by the user, these angles are obtained from the 
spline fits through the blade surface coordinates. Since spline curves pass exactly 
through the fitted points (see refs. 7 and 8), these points must be given accurately or 
else a "wavy" curve will result. And the angles on a "wavy” curve are greatly in error. 
Since the forces computed in this program depend heavily on these angles, it is impor- 
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tant that the curve fits be accurate. The user should always check the spline curve out- 
put to the program. If a particular curve is not smooth (i. e. , its angles do not increase 
or decrease smoothly, or the curvature changes sign radically from point to point), one 
of two changes should be made. Either the input points should be adjusted so that the 
curve fit improves, or input angles should be supplied in that particular region by the 
user. In the numerical example, input angles could be given for the relatively poor 
spline angles (47.26 and 40.37) at the third input points on surfaces 1 and 3 of example 
tandem blade number 2 (see table IV). 

TABLE IV. - SELECTED OUTPUT FOR TANDEM BLADE NO. 2 OF NUMERICAL EXAMPLES 


TANDEM OCA BLADE C2/C1=1.0 PH2/PH1=2.0 G/Cl=.l L/C1=.2 F=1.4 


GAM 

R 

PTZ 

TT 2 



1.400 

1116.48 

2080.00 

520.00 



VR I 

VRO 

BETAI 

BETAC 

ALPHMX 


51.2910C 

249.63800 

43.250 

-1.000 

0. 


KPV 

KSURF 

NSURF1 

NSURF2 

NSURF3 

NSURF4 

2 

4 

27 

36 

28 

36 


SLCTION 

SURFACE - FRONT 

BLADE 





XI 

Y 1 

VI 

ANG1 




c. 

0. 

389.0000 

-0. 




C .00040 0 

0.0C1150 

389.7000 

-0. 




C.001CCC 

0.0C1920 

3 90.5000 

-0. 




C.002C0C 

0.003240 

392.3000 

-0. 




C. 005000 

0.007050 

397.1000 

-0. 




0.015164 

0.019214 

411.2000 

-0. 




C. 025274 

0.029916 

422.7000 

-0. 




C. 035384 

0.039514 

430.0000 

—0. 




SURFACE 

1 - SPLINE FIT 

OF X-Y 



SPL INE 

INPUT 

X 

Y 

SLOPE 

SEC. DERIV. 

CURVATURE 

ANGLE 

ANGLE 

C. 

0. 

4.06782 

-7156.91138 

-97. 36791 

76.1887 

-0. 

C.00040C 

0.001150 

1.92075 

-3578.45569 

-352.40246 

62.4971 

-0. 

C. 001000 

0. OC 1 920 

1.08205 

782.79,342 

244. 74290 

47.2567 

-0. 

C.00200C 

0.0C3240 

1.40451 

-137.86(566 

-26.89928 

54*5494 

-0. 

C .005000 

0.0C7050 

1.20778 

6.71066 

1. 74060 

50*3764 

-0. 

0.015164 

C. 019214 

1.14066 

-19.91-772 

-5.70605 

48.7554 

-0. 

0.025274 

0.029916 

0.99503 

-8.89076 

-3. 16687 

44.8573 

-0. 

0.035384 

0.039514 

0.90295 

-9.32564 

-3.81286 

42.0804 

-c. 

0.045493 

0.048186 

0. 814 79 

-8. 116 54 

-3.78171 

- 39. 1727 

-0. 

0.055603 

0.056028 

0. 73846 

-6.98326 

-3.63522 

36.4443 

-0. 

C. 065712 

0.063142 

0.66936 

-6.68513 

-3.83650 

33. 7969 

-c. 

C. 075822 

0. 069577 

0.60475 

-6.09907 

-3.82139 

31.1633 

-0. 

C. 085932 

0.075391 

0. 54656 

-5.41140 

-3.65622 

28.6593 

-0. 

0.096041 

0.08C644 

0.49314 

-5.15797 

-3.72113 

26.2497 

-c. 

C. 106150 

0. 065368 

0.44171 

-5.01593 

-3.83919 

23.8317 

-0. 

C.11626C 

0.089584 

0.39297 

-4.62750 

-3.73078 

21.4532 

-0. 

0.12637C 

0.093328 

0.34844 

-4.18121 

-3. 52092 

19.2103 

-0. 

C.13648C 

0.096622 

0.30171 

-5.06328 

-4.44302 

16. 7850 

-0. 
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TABLE TV. - Continued. SELECTED OUTPUT FOR TANDEM BLADE NO. 2 OF NUMERICAL EXAMPLES 


C .143 130 

0.098534 

0.27598 

-2.67535 

-2.39642 

15.4283 

-0. 

0.14977C 

C. 10C268 

0. 24036 

-3.05169 

-7.40112 

13.5154 

-c. 

C .1520GC 

C. 10C 81 C 

0.25740 

23.33112 

21. 19065 

14.4346 

-0. 

C.1531CC 

C. 101078 

0.20328 

-121.73237 

-114.55878 

11 .4904 

-0. 

C.1540CC 

0. 101270 

0.28822 

310.50212 

275.47210 

16. 0783 

-0. 

C.15500C 

0. 1C1500 

-0. 041 70 

-970.35752 

-967.83167 

-2.3880 

-0. 

0.15580C 

0. 1 Cl 300 

-0.27845 

378.49150 

338.38385 

-15.5598 

-0. 

C .156200 

C. 1CC930 

-2. 29381 

-10455 .30090 

-667. 28804 

-66.4449 

-0. 

0.156420 

C. 100088 

-5. 74407 

-20910.60181 

-105.50107 

-80.1242 

-0. 


SURFACE 2 

- SPLINE FIT 

OF X-Y 



SPL INE 

INPUT 

X 

Y 

SLOPE 

SEC. DERIV. 

CURVATURE 

ANGLE 

ANGLE 

C. 

0. 

-4.04298 

7907.87189 

109.46437 

-76.1072 

-0. 

0 .000400 

-0.001090 

-1.67062 

3953.93594 

535.68018 

-59.0961 

-0. 

C.00100C 

-O.OClfcOO 

-0.39495 

298.29909 

240.00753 

-21.5514 

-0. 

C .002000 

-0.001800 

0.04074 

573.08138 

571.65733 

2.3331 

-C. 

C .003000 

— C.0C1470 

0.62197 

589.37535 

360.86779 

31.8805 

-0. 

0 .00400 C 

-0. 0CC670 

0.86137 

-110.58288 

-48.09860 

40.7406 

-0. 

C. 005000 

0. 0CC150 

0. 79255 

-27.04368 

-13.01748 

38.3988 

-0. 

C .006000 

0.0CC540 

0. 79841 

38.75761 

18.49696 

38.6043 

-0. 

C.0C800C 

0.0C2580 

0.82442 

-12.75106 

-5.85749 

39.5028 

-0. 

C. 015164 

0.008266 

0.77791 

-0.23308 

-0. 11461 

37.8797 

-0. 

C .025274 

0. 016086 

0. 76583 

-2. 15563 

-1.07873 

37.4461 

-0. 

0.035384 

0.023694 

0. 73679 

-3.58913 

-1. 87282 

36.2826 

-0. 

0.045493 

0.031025 

0.72014 

0.29372 

0. 15695 

35.7591 

-c. 

C. 055603 

0. 038250 

0. 70216 

-3.85091 

-2. 11088 

35.0749 

-c. 

C .065713 

0.045199 

0.67717 

-1.09174 

-0.61978 

34.1047 

-0. 

0.075822 

0.051956 

0.65642 

-3.01410 

-1. 76096 

33.2816 

-0. 

C .080 87 7 

0.055250 

0.64968 

0.34737 

0.20484 

33.0109 

-0. 

0.085532 

0.058522 

0.64161 

-3.54143 

-2.11147 

32.6844 

-0. 

C .090987 

0.061731 

0.63019 

-0.97438 

-0.59002 

32.2188 

-0. 

C. 096041 

0. 064898 

0.62198 

-2.27666 

-1.39396 

31.8806 

-0. 

0.101100 

0.068021 

0.613 75 

-0.97385 

-0.60288 

31.5397 

-0. 

C . 106 15 C 

0. C 71 102 

0.60525 

-2.39509 

-1.49966 

31.1043 

-0. 

0.111210 

0.074120 

0.58489 

-5.65009 

-3. 63394 

30.2231 

-0. 

C. 114000 

0. 075750 

0.59078 

9.87256 

6. 30086 

30.5739 

-0. 

0.116260 

0.077095 

0.59268 

-8. 19776 

-5. 21896 

30.6542 

-0. 

C. 121320 

0.080028 

0.57432 

0.94288 

0.61483 

29.8697 

-0. 

0.126370 

0.082918 

0.56581 

-4. 31413 

-2.84422 

29.5016 

-0. 

C .131430 

0. 085744 

0. 55479 

-0.04071 

-0.02722 

29.0211 

-c. 

C.13648C 

C. 088549 

0. 55686 

0.85919 

0. 57298 

29.1116 

-0. 

0. 143130 

0.092210 

0.53501 

-7.43088 

-5.09402 

28.1471 

-0. 

C. 149770 

0.055752 

0. 55496 

13.44039 

8.98485 

29.0284 

-0. 

0.153000 

0.097480 

0.47333 

-63.98429 

-47. 24755 

25.3296 

-0. 

0.155000 

0.098440 

0. 55733 

147.97984 

98. 62603 

29.1320 

-0. 

C.15580C 

0.098650 

0.36366 

-632. 15340 

-524. 70058 

19.9841 

-c. 

C .156200 

0.099240 

2.32413 

10434. 53406 

644. 23296 

66.7193 

-0. 

0. 15642C 

0. 100C88 

5.76754 

20869.068 12 

104.04876 

80.1636 

-0. 


surface 2 - 

SPLINE FI T 

OF X-Y 



SPL INE 

INPUT 

X 

Y 

SLOPE 

SEC. OERIV. 

CURV ATURE 

AN GLE 

ANGLE 

0. 136480 

C. 057651 

4. 56656 

-10132.052 12 

-99.17828 

77.6482 

-0. 

0.136800 

0.05868C 

2.13486 

-5066.02606 

-386. 66180 

64.9010 

-0. 

C.1374CC 

0. 059400 

0. 85009 

783.42208 

346.49924 

40.3674 

-0. 

C. 138000 

0.060000 

1. 06480 

-67.71216 

-21. 72389 

46.7975 

-0. 

C.13980C 

0. 061 838 

0.99468 

-10.19747 

-3.63424 

44.8472 

-0. 

C. 146450 

0.068170 

0. 90109 

-17.95037 

-7.35954 

42.0216 

-0. 

0.15310C 

0. 073822 

0.80728 

-10.26168 

-4.83414 

38.9133 

-0. 

0.161460 

C.08C198 

0.71637 

-11.48868 

-6. 17211 

35.6165 

-c. 
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TABLE IV. - Continued. SELECTED OUTPUT FOR TANDEM BLADE NO. 2 OF NUMERICAL EXAMPLES 


SURFACE A - 

SPLINE FIT 

OF X-Y 



SPL INE 

INPUT 

X 

Y 

SLOPE 

SEC. OERIV. 

CURV ATURE 

ANGLE 

ANGLE 

0.136480 

0. 057651 

-4.40039 

10245.14673 

111.49134 

-77.1568 

-0. 

C.1368CC 

0.056680 

-1.94155 

5122.57336 

491.76962 

-62.7492 

-0. 

0.137400 

0.056120 

-0.45367 

-162.95498 

-123.06826 

-24.4023 

-0. 

C.13800C 

0. 055870 

-0.29377 

695.94306 

614.67782 

-16.3714 

-0. 

C.13900C 

C.05590C 

0.32957 

550.75204 

471.82649 

18.2409 

‘0. 

C. 140000 

0.056430 

0. 65547 

101.04542 

59.11178 

33. 2438 

-0. 


SURFACE 1 - 

INTEGRATION OF 

(AVERAGE PRESSURE*TAN( AVERAGE 

ANGLED TO GET FX1 


X 

PRE SSURE 

VELCCI TY 

I NT EGRAN01 

INT EGRAND2 

SUM 

C . 

1908.95782 

389.00000 

0 . 

5063.37781 

C. 

C. 000 400 

1908. 36055 

389.70000 

5061.79358 

2712.99512 

2.02503 

C.001000 

1907.67673 

390.50000 

2712.02298 

2347.65125 

3.65254 

C .002000 

1906.13377 

392.30000 

2345.75244 

2480. 79465 

5.99924 

C .005000 

1901.98886 

397.10000 

2475.40012 

2232.29367 

13.43353 

C. 015164 

1889.56056 

411.20000 

2217.70703 

2012.76736 

36.04844 

0.025274 

1879.14795 

422. 70000 

2001.67581 

1781.30598 

56.34145 

C .035384 

1872.41116 

430.00000 

1774.91997 

1606.35677 

74.31817 

C. 045493 

1865.24011 

433.40000 

1603.63629 

1450.37564 

SC. 54308 

0 .055603 

1870. 54842 

432.00000 

1451.39078 

1315.64914 

105.21151 

0.065713 

1873.33980 

429. 00000 

1317.61246 

1192.53510 

1 18.52264 

C. 075822 

1877.95544 

424.00000 

1195.47333 

1080 . 36 1 82 

130. 59283 

0.085932 

1883. 70537 

417.70000 

1083.67268 

978.69334 

141.53204 

0.09604 1 

1890.27571 

410.40000 

982.10702 

883.08479 

15 1.44290 

0.106150 

1898.58884 

401.00000 

886 .96844 

791.95667 

160.38964 

0. 116260 

1907.67673 

390.50000 

795.74750 

706.87220 

168.41548 

C .1263 70 

1918.30307 

377.90000 

710.80969 

623.28059 

175.58186 

0.136480 

1928. 70177 

365.20000 

626.65925 

557.00500 

181.90031 

0.143130 

1937.84424 

353.70000 

559.64532 

500.14394 

185.61317 

0 . 149 770 

1947.48326 

341.20000 

502 .63171 

484.65864 

188.94238 

0 .152000 

1950. 72658 

336.90000 

485.46579 

449.01623 

190.02407 

C .153 IOC 

1952.14780 

335.00000 

449.34336 

478.92900 

190.51817 

C . 15400 C 

1953. 63586 

333 .00000 

479.29408 

234.51831 

190.94937 

C. 155000 

1955. 11591 

331.00000 

234.69597 

-308.74772 

191.18397 

0 . 15 5 80 C 

1956. 36755 

329.30000 

-308.94538 

-1700.78606 

ISC. 93689 

0.156200 

1956.95448 

328.50000 

-1701.29631 

-6516.48517 

190.25648 

0 . 156 42 C 

1957. 32066 

328.00000 

-6517.70447 

0 . 

188.82271 


SURFACE 1 - 

INTEGRATION OF 

(-P) TC GET 

FY1 



X 

PRESSURE 

VELOCI TY 

I NT EGRAN01 

INT £ GRAN 02 

SUM 

C. 

1908.95782 

389.00000 

-1908 .95782 

-1908.95782 

0 . 

C. 000400 

1908.36055 

389.70000 

-1908. 36055 

-1908.36055 

-C. 76346 

0 .001000 

1907.67673 

390.50000 

-i907. 67673 

-1907.67673 

- 1.9C827 

C. 002000 

1906.13377 

392.30000 

-1906.13377 

-1906. 13377 

-2.81518 

0 .005000 

1901. 98886 

397.10000 

-1901.98886 

-1901.98386 

-9.52736 

C. 015164 

1889. 56056 

411.20000 

-1889.56056 

-1889.56056 

- 28.79602 

0.025274 

1879.14795 

422.70000 

-1879.14795 

-1879.14795 

-47.84684 

C .035384 

1872.41116 

430.00000 

-1872.41116 

-1872.41116 

-66.81097 


SURFACE 2 - 
X 

INTEGRATION GF 
PRE SSURE 

(-AVERAGE 
VELOCI TY 

PRESSURES AN( AVERAGE 
I NT E GRAND1 

ANGLE)) TO GET FX 2 
INT EGRAND2 

SUM 

C. 

1986.24275 

286.00000 

0 . 

4819.36914 

0 . 

0.000400 

1986. 37170 

285.80000 

4819.68201 

1685.97951 

1.92781 

C.OOIOOC 

1986. 56491 

285.50000 

1686.14349 

336.32727 

2.93945 

C. 002000 

1986. 88658 

285.00000 

336.38173 

-611.50432 

3.27580 
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TABLE IV. - Continued. SELECTED OUTPUT FOR TANDEM BLADE NO. 2 OF NUMERICAL EXAMPLES 


C. 003000 

1987.14349 

284.60000 

-611.58339 

-1460.26328 

2.66426 

0.004000 

1987.40013 

284.20000 

-1460.45187 

-1642.34860 

1.2C390 

0.00500C 

1987. 59236 

283.90000 

-1642. 50746 

-1581.08849 

-0.43853 

C.0060QC 

1987.91229 

283.40000 

-1581.34297 

-1612.85777 

-2.C1974 

c.ooaooo 

1988.23177 

282.90000 

-1613.11697 

-1592.37494 

-5.24572 

0.01516 4 

1989.44083 

281 .00000 

-1593.34328 

-1535.55826 

- 16.65696 

0.025274 

1991.27148 

278.10000 

-1536.97127 

-1495.86957 

-32.18860 

C. 035384 

1994.20064 

273.40000 

-1498.06998 

-1452.64209 

-47.32296 

0.045493 

1997. 38667 

268.20000 

-1454.96291 

-1420.35960 

-62.01945 

0.055603 

2000.27650 

263.40000 

-1422.41458 

-1379.37242 

-76.38967 

C. 06571 3 

20C3. 17566 

258.50000 

-1381.37166 

-1335.60754 

-90.24523 

C. 075822 

2006. 36848 

253.00000 

-1337.73634 

-1310.24495 

- 103.85765 

C. 080877 

20C8. 30965 

249.60000 

-1311.51262 

-1296.63350 

- 110.48414 

C. 085932 

2010. 00174 

246.60000 

-1297.72597 

-1278. 12712 

- 117.04138 

C. 090987 

2011. 84050 

243.30000 

-1279.29636 

-1259.56548 

- 123.50527 

0.096041 

2013. 76498 

239.80000 

-1260.77036 

-1244.21977 

-129.87416 

0.101100 

2015.33940 

236.90000 

-1245.19255 

-1228.33456 

- 136.17113 

C. 106150 

2016. 521 59 

234.70000 

-1229.055 10 

-1199.88022 

- 142.37604 

C. 111210 

2017. 74619 

232.40000 

-1200.60889 

-1186.10077 

- 146.44927 

C. 114000 

2017.95802 

232.00000 

-1186.22530 

-1194.08492 

-151.75867 

C. 116260 

2017.85213 

232.20000 

-1194.02226 

-1177.33980 

- 154.45723 

0.121320 

2016. 73538 

234.30000 

-1176.68823 

-1149.65436 

- 160.41292 

0.12637C 

2013. 71034 

23 9.90000 

-1147.92992 

-1128.25554 

-166.21432 

0.131430 

2009.55240 

247.40000 

-1125.92589 

-1116.95712 

- 171.91739 

C .13648 C 

2003. 17566 

258.50000 

-1113.41278 

-1093.49748 

-177.54908 

C. 143130 

1991. 33432 

278.00000 

-1087.03351 

-1085.15915 

- 184.79934 

C .149770 

1977.00883 

300.00000 

-1077.35260 

-1015. 12959 

-191.97888 

C. 153000 

1968. 13889 

312.90000 

-1010.57516 

-1012.82432 

-195.25C39 

0.155000 

1962.32207 

321.10000 

-1009.83092 

-896.68433 

-157.27305 

C.15580C 

1959.65224 

324.80000 

-895.46435 

-1850.02377 

- 197.98991 

0.156200 

1958.19746 

326.80000 

-1848.65039 

-6586.04169 

- 198.72964 

0.15642C 

1957. 32066 

328.00000 

-6583.09271 

0 . 

-200. 17825 


SURFACE 2 - 
X 

integration of 

PRESSURE 

IP) TO GET FY2 
VELCCI TY 

INTEGRAND1 

INT EGRAN02 

SUM 

0 . 

1986.24275 

286.00000 

1986.24275 

1986.24275 

0 . 

C.C004QO 

1986. 371 70 

285.80000 

1986.37170 

1986.37170 

C. 79452 

0.001000 

1986. 56491 

285.50000 

1986.56491 

1986.56491 

1.98640 

0.002000 

1586. 88658 

285.00000 

1986.88658 

1986.88658 

3.97313 

0 .003000 

1987.14349 

284.60000 

1987.14349 

1987. 14349 

5.96C14 

0 .004000 

1987.40013 

284.20000 

1987.40013 

1987.40013 

7.94742 

0.005000 

1987. 59236 

283.90000 

1987.59236 

1987. 59236 

9.93491 


SURFACE 3 - 

INTEGRATION OF 

(AVERAGE PRESSURES AN< AVERAGE 

ANGLE)) TO GET FX 3 


X 

PRESSURE 

VELOCl TY 

INTEGRAND1 

INT EGRAND2 

SUM 

C.13648C 

2025.15070 

218.00000 

0 . 

5974.31323 

0 . 

C. 136800 

2024. 24997 

219.30000 

5971.65601 

2650.88300 

1.91136 

C.1374Q0 

2022.63115 

223.00000 

2648.76303 

1924.94296 

3.50125 

C .138000 

2020. 98991 

226.20000 

1923.38100 

2079.85095 

4. 65575 

C. 139800 

2016.36099 

235.00000 

2075.08719 

1909.07281 

8.39519 

C.14645C 

2000.44951 

263.11000 

1894.00795 

1706.58331 

21.04044 

0 . 153 IOC 

1985.33784 

287.40000 

1693.69154 

1510.50345 

32.34635 

C. 161460 

1968. 76892 
/ 

312.00000 

1497.8973 2 
1 

1 304. 13586 

44.92146 
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TABLE IV. - Continued. SELECTED OUTPUT FOR TANDEM BLADE NO. 2 OF NUMERICAL EXAMPLES 


SURFACE 2 - 

INTEGRATION OF 

(-P) TG GET 

F Y3 



X 

PRESSURE 

VELOCI TY 

I NT EGRAND 1 

INT EGRAND2 

SUM 

C. 136480 

2025.15C70 

218.00000 

-2025.15070 

-2025.15070 

0 . 

C. 136800 

2024.24997 

219.80000 

-2 024.24997 

-2024.24997 

-0.64791 

C.13740C 

2022.63115 

2 23.00000 

-2022.63115 

-2022.63115 

- 1.86197 

C. 138000 

2020.98991 

226.20000 

-2020.98991 

-2020.98991 

-3.07506 

C. 139800 

2016.36099 

235.00000 

-2016.36099 

-2016.36099 

-6.70867 

C .146450 

2000.44951 

263.1 1000 

-2000.44951 

-2000.44951 

-20.06457 

C. 153100 

1985.33784 

287.40000 

-1985.33784 

-1985.33784 

-33.31731 

C. 161460 

1966. 76892 

312.00000 

-1968.76892 

-1968.76892 

-49.84547 


SURFACE A - INTEGRATION OF (-AVERAGE PRE SSURE*T AN ( AV E R AGE ANGLE) I TO GET FX4 


X 

PRESSURE 

VELOCI TY 

INT EGRAND1 

INT EGRAND2 

SUM 

C .136 480 

1823.69095 

480.00000 

0 . 

5003.21576 

0 . 

0.136800 

1830.77676 

473.00000 

5022.65540 

1741.94926 

1.6C414 

0.137400 

1845.67378 

458.00000 

1756.12350 

685.91782 

2.65356 

0.138000 

1860.64937 

442.50000 

691.48328 

-30.35918 

3.06678 

C. 139000 

1883.07016 

418.40000 

-30.72501 

-907.97530 

3.03624 

C. 140000 

1901.032 82 

398.20000 

-916.6365 1 

-1267.87556 

2.12394 

0 .141000 

1917.22090 

379.20000 

-1278.67204 

-1256.66789 

C. 85066 

C. 143000 

1945.73315 

343.50000 

-1275.35663 

-1205.80515 

- 1.68136 


SURFACE A - INTEGRATION OF (P) TO GET FY4 


X 

PRESSURE 

VELOCITY 

I NT EGRAND1 

INT EGRAND2 

SUM 

C. 136480 

1823.69095 

480.00000 

1823.69095 

1823.69095 

0. 

C. 136800 

1830.77676 

473.00000 

1830.77676 

1830.77676 

C. 58472 

C. 137400 

1845.67378 

458.00000 

1845.67378 

1845.67378 

1.68765 

C. 138000 

1860.64937 

442.50000 

1860.64937 

1860.64937 

2.79955 

C. 139000 

1883.07016 

418.40000 

1883.07016 

1883.07016 

4.67141 

0 .14000 C 

1901.032 82 

398.20000 

1901.03282 

1901.03282 

6.56346 

C. 141000 

1917.22090 

379.20000 

1917.22090 

1917.22090 

8.47259 

C.14300C 

1945. 73315 

343.50000 

1945.73315 

1945.73315 

12.33554 


ANGLE BETWEEN HER IOIONAL- TANGENT IAL AXES ANO LIFT-DRAG AXES = 25.0593 


SUCTION SURFACE FORCES - FRONT BLADE 

X Y MERIDIONAL TANGENTIAL DRAG LIFT 

FORCE FORCE FORCE FORCE FORCE FORCE 

188.8227 -296.5217 188.8227 -296.5217 45.4553 -348.5871 

PRESSURE SURFACE FORCES - FRONT 8LACE 

X Y MERIDIONAL TANGENTIAL DRAG LIFT 

FORCE FORCE FORCE FORCE FORCE FORCE 

-200.1762 312.9505 -200.1782 312.9505 -48.7834 368.2791 


23 



TABLE IV. - Concluded. SELECTED OUTPUT FOR TANDEM BLADE NO. 2 OF NUMERICAL EXAMPLES 


SUCTION SURFACE FORCES - REAR BLADE 


X 

Y 

MERIDIONAL 

TANGENTIAL 

DRAG 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

84.2569 

-352.4448 

84.2969 

-352.4448 

-72.9131 

PRESSURE SURFACE 
X 

FORCES - REAR 
Y 

BLADE 

MERIDIONAL 

TANGENTIAL 

DRAG 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

-87. Cf 51 

365. 9404 

-87. 0851 

365.9404 

76.1085 


SUM OF SUCTION 

SURFACE AND 

PRESSURE SURFACE 

FORCES - FRONT 

EL ADE 

X 

Y 

MERIDIONAL 

TANGENTI AL 

DRAG 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

-11.3555 

16.4288 

-11.3555 

16.4238 

-3.3281 

TOTAL FORCE ON FRONT BLADE = 19.9713 

ANGLE WITH RESPECT TO MERIDIONAL AXIS = 

124.6521 



SUM OF SUCTION SURFACE AND PRESSURE SURFACE 
X Y MERIDIONAL 

FORCE FORCE FORCE 

-2. 7682 13.4956 -2.7882 


FORCES - REAR BLADE 
TANGENT I AL 
FORCE 
13.4956 


DRAG 

FORCE 

3.1904 


TOTAU FORCE ON REAR BLADE = 13.7806 

ANGLE HUH RESPECT TO MERIDIONAL AXIS = 101.6732 


RATIO OF HEAR BLADE 
X 

FORCES 

0.2455 


FORCES TO FRONT BLADE FORCES 
Y MERIDIONAL 

FORCES FORCES 

0.8215 0.2455 


TANGENT IAL 
FORCES 
0.8215 


DRAG 

FORCES 

-0.9586 


RATIO OF TOTAL FORCE ON REAR BLADE TO TOTAL FORCE ON FRONT BLADE = 


SUM OF ALL FRONT 

ANO REAR 8LADE FORCES 



X 

Y 

MERIDIONAL 

TANGENT IAL 

DRAG 

FORCE 

FORCE 

FORCE 

FORCE 

FORCE 

- 14. 1438 

29.9244 

-14.1438 

29.9244 

-0.1377 

TOTAL FORCE 

ON ENTIRE 

TANOEM BLADE = 

33.0986 


ANGLE WITH 

RESPECT TO 

MERIDIONAL AXIS = 

115.2977 



L IF T 
FORCE 
-354.5736 


LIFT 

FORCE 

368.3798 


LIFT 

FORCE 

19.6921 


L IF T 
FORCE 
13.4062 


L IFT 
FORCES 
0 .6808 

C.650C 


LIFT 

FORCE 

33.0983 
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OUTPUT 


Table IV contains a portion of the output generated for tandem blade number 2 of the 
numerical examples. Most parts of the output have been abbreviated due to their length. 

The first portion of output contains a copy of all the input to the program. The input 
is shortened here because it was contained in its entirety in table II. 

Following the input are the spline curve fits of the X-Y input coordinates of each of 
the blade surfaces. The X and Y coordinates are listed, followed by the slope, second 
derivative, and curvature at each point on the surface. From the slopes, angles are 
computed along the spline curves. These are also printed, along with any input angles 
which were given by the user. The user should always check the angles and second de- 
rivatives of the spline curve to see if a good fit has been obtained. In regions where the 
spline fit is poor (fluctuating angles or strongly fluctuating second derivatives), the input 
data should be altered for a subsequent run on the program. Either input data points 
should be shifted slightly, or points added or subtracted in the critical regions, or angles 
given as input to be used instead of the spline angles. 

The next portion of output contains the summation of components of pressure force 
along the blade surfaces (see eqs. (2) to (9)) to obtain force contributions in the X and 
Y directions. A force in the X direction and a force in the Y direction are computed 
on each of the blade surfaces. The X coordinates of the surface points are listed, to- 
gether with the static pressure and velocity at the points. The columns labeled "Inte- 
grand 1" and "Integrand 2" contain the quantities 


±P^ tan 


+ a 


i-1 


2 


and 


±P i tan 


0!j + O’ 


i+1 


2 


for the F x integrations, and the quantities iP^^ for the Fy integrations. The column 
"Sum" contains the cumulative summation of force at each point along the surface, so 
that the final value of this column is the total force on a surface in the X or Y direc- 
tion. The differences in the elements in this column indicate the contribution to force 
from each element of the surface, and help to locate errors in the input when they exist. 
The final portion of output lists all the resultant forces on blade surfaces and whole 

blades. At the beginning of this section, „„„ (the angle between the M-0 axis and the 

iiiGcin 
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lift-drag axis) is listed. This is followed by lists of forces in the X and Y directions, 
the M and Q directions, and the lift and drag directions. For a tandem blade, the lift- 
drag direction is perpendicular and parallel to the mean velocity vector across the over - 
all blade, not across each segment individually. Thus, drag forces listed for the indi- 
vidual segments are not really in the true drag directions for those segments. 

Components of force are given first for each of the blade surfaces , and are then 
combined to give forces on each of the individual blade segments . Total forces are also 
given on the blade segments and are located with respect to the M or Z axis. For 
tandem blades , ratios are given for all rear segment forces to front segment forces. 
Finally, for tandem blades, forces are summed from all four surfaces to give total 
forces bn the overall tandem blade section. 

COMPLETE PROGRAM LISTING 


SIBJOB 

tIBFTC LIFTOR 

C 

c liftdr computes forces due to pressure on either a single 

C OR TANDEM BLADE. GIVEN SURFACE GEOMETRY AND SURFACE PRESSURE 
C OR VELOCITY DISTRIBUTIONS 

f. 

DIMENSION Xl!LOO),X21 100) » X3 4 100) , X4!IQQ)»Y1 1100), Y2! 133). 

1Y3! 100) .Y4I 100) ,P1! 100) ,P2< 103) .P31100) , P4( 130 ) , VI 1 100 ) , V 21 100), 
2V3( 100) ,V4I 100) ,ANG1 ( 100) . ANG2 1100 ) , ANG3 C 100 ) , ANG4 ( 100 ) , 

3ANGS1I 100).ANGS2(100), ANG S3 ( 100) .ANGS4I100) , SL 1 1 00 ) , SDI 1 3 3 ) . 

4CU! 100) ,Z ( IOC. 2) .SUM ( 100) 

C 

C READ AND PRINT INPUT 
C 

10 WRITE! 6.1000) 

READ IS. 1040) 

WRITE! 6. 1040 ) 

REAO 1 5,1030) GAM. R.PTZ.TTZ 
WRITE! 6. 1050) GAM, R.PTZ.TTZ 
READ 15.1030) VRI , VRO.BETAI , BE TAO, ALPHMX 
WRITE! 6,1060) VRI, VRO.BETAI , BE TAG. AL PHMX 
REAO 15.1010) KPV.KSURF .NSURE 1 .NSURF2.NSURF3, NSURF4 
WRITE! 6,1070 ) KPV, XSURF , NSURF 1 .NSURF2 , NSURF3 , NSURF 4 
IF 1 1XPV.NE. 1. AND.KPV.NE. 2). OR. 1 XSURF.NE.2. AND.KSURF.NE.4) ) GO TO 
1320 

IF 1KSURF.E0.2) WR I TE ( 6, 1080) 

IF (KSURF.E0.4) WR I TE 1 6. I 1 00 ) 

IF (KPV.EO.l) WRIIE16.1 140 ) 

IF 1KPV.E0.2) WRITE16.il 50) 

READ 1 5.1020) 1 XU I ) ,Y1 II) , VI ( I ) , ANG1 1 1 ) . 1 = 1 .NSURF l) 

WRITE! 6,1220 ) ( X 1 1 I ) , Y 1 ( I ) , VI { I ) • ANG1 i I ) .1 = 1 .NSURF 1) 

IF (KSURF.E0.2) WR I TE 1 6, 1 090) 
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n n n r> r> 


IF (KSURF.E0.4) WR I TE 1 6, 1 1 101 
IF (KPV.E0.1) MR I TE ( 6 • 1 160 ) 

IF (KPV.E0.21 WRITE(6, 11701 

READ ( 5.10201 ( X2( I 1 . Y2( 1 1 . V2( 1 1 ,ANG2 ( 1 1 .1 = 1 .NSURF2) 

WRITE! 6.1220 1 ( X2 ! I ) « Y2 ( I 1 « V2 ( 1 ) . ANG2 ( I 1.1 = 1 .NSURF2) 

IF (KSURF.E0.2I GO TO 20 
MRI TE( 6.1120) 

IF (KPV.E0.11 WRITE(6, 11801 
IF (KPV.E0.21 HRITEI6. 11901 

READ ( 5*10201 ( X3( I 1 . Y3 ( I 1 . V3 ( 1 1 . ANG3 ( I 1 , 1= 1 . NSURF 3 1 
WRITE! 6.1220 1 ( X3( I 1 . Y 3 ( I 1 , V3 ( 1 1 .ANG3 ( I 1 ,I = 1.NSURF31 
WRITE! 6. I 130 1 

IF (KPV.EO.l) WRITE(6. 12001 
IF (KPV.E0.21 WRITE (6.1210) 

READ ( 5.1020) ( X4( 1 1 , Y4( I 1 , V4( 1 1 , ANG4 ( I 1 , 1 = 1 , NSURF4 1 
WRITE! 6.1220) ( X4( I) . Y4( I 1 , V4 ( I 1 , ANG4 ( I 1 ,I=1,NSJRF4) 

20 8 ET A 1= 8ETAI /57. 295780 
8ETA0= BETAO/57. 295780 
ALPHMX= ALPHMX/57. 295780 
C 

C FIT SPLINE CURVES TO X AND Y. COMPUTE SURFACE ANGLES AND 
CURVATURES, AND PRINT RESULTS 

CALL SPLINE! XI . Yl , NSURF 1 , SL .SO ,CU, ANGS1 1 
J = 1 

WRI TE( 6. 1000 1 
WRITE! 6,1230) J 

WRITE! 6.1240) ( XI ( I ) . Y 1 ( I 1 . SLI 1 1 , SD ( I 1 ,CU( 1 1 , ANGS 1 ( Il.ANGK I 1 
1 1=1. NSURF 1) 

C AL L SPL I NE ( X2 * Y2 , NSURF2 . SL . SO ,CU. ANGS2 1 

J = 2 

WRITE! 6.1230) J 

WRITE! 6. 1240) ( X2( I 1 ,Y2( I ) . SLU ) .SD! I ) ,CU( 1 1 , ANGS2 ( I ) . ANG2 ( I 1 
II =1.NSURF2) 

IF (KSURF.E0.2) GO TO 30 

CALL SPLINE! X3, Y3. NSURF 3, SL. SO, CU.ANGS3) 

J= 3 

WRITE! 6.1230) J 

WRITE! 6. 1240) ( X3( 1 1 , Y3 ( 1 1 . SLi 1 1 • SD ( I 1 .CU( I ) , ANGS3 ( 1 1 « ANG3! I) 
1 1 =1, NSURF 3 1 

CALL SPLINE! X4. Y4.NSURF4, SL .SD.CU. ANGS4) 

J = 4 

WRITE! 6.1230) J 

WRITE! 6.12401 ( X4( 1 1 . Y4( I 1 , SL! 1 1 ,SD( I I ,CU! 1 1 , ANGS4 ( I), ANG4! I) 
1 1 = 1. NSURF 4 1 

COMPUTE VELOCITY AND PRESSURE DISTRIBUTIONS ALONG SURFACES 

30 R HTZ = PTZ/R/TTZ 

T 1= (GAM- 1. 1 /GAM*RHTZ/PTZ/2. 

T 2= l./Tl 
T 3= GAM/! GAM-1. I 
T 4= 1./T3 

IF (KPV.EO. 1 ) GO TO 50 
DO 40 1=1. NS UR FI 

40 P1(I)= PTZ*( 1.-T1*V1U)**2)**T3 
GO TO 70 

50 DO 60 I =1 .NSURF 1 
PI! I)= VI! I) 



u u o a 


60 V 1C I »= SQRT(T2*(l.-(Pl(I)/PTZI**T41) 

GO TO 90 

70 00 80 l =1 .NSURF 2 
80 P ?( I ) = PTZ*( l.-n*V2(I >**21 **T3 
GO TO 110 

90 00 LOO 1= l.N SURF 2 
p n i >= V2( i ) 

100 V 2( I ) = SORT! T2*(l.-(P2(II/PTZ)**T4)I 
110 IF IKSUkF.eO.2l GO TO 190 
IF (KPV.EO.l » GO TO 130 
00 120 I=1,NSURF3 

120 P 31 11= PTZ*( l.-Tl*V3( I )**2)**13 
GO TO 150 

130 00 1£0 I = 1. N SURF 3 
P 31 l ) = V3( I) 

140 V 31 11= S0RTII2*C1.— CP3II WPTZI 4*14) ) 

GO TO 170 

150 00 160 I = i,NSURFA 

160 P 4( I ) = PTZ*( l.-Tl*V4( I )**2)**T3 

GO TO 190 

170 00 180 I = 1 * N SURF4 
P -4C I ) = V4( I ) 

180 V 4( I ) = SORT ( T2* (1.-4P44 I ) / PTZ1 **T4) ) 

INTFGRATF PRESSURE DISTRIBUTIONS TO OBTAIN COMPONENTS OF X AND Y 
FORCES. AND PRINT RESULTS 

190 NM= NSURF 1-1 

IF ( ANGH 1). £0.0.1 ANGHl)=ANGSnn 
IF ( ANGH 2I.FO.O.I ANG 1 ( 2 » =ANS Sl( 2 ) 

ANG 1(1)— ANGllll/57. 295780 
ANG1I 21= ANGK2I/57. 295780 
Z Cl. 11= 0. 

Z (1.21= Pit 1I*TAN( ( ANG 1 ( 1 ) + ANG 1(2) 1/2.) 

00 200 1=2. NM 

IF ( ANG1( I + U.EO.O.I ANG1 ( mi = ANG SI (HI ) 

ANGH IH) = ANG1( I*-ll/57. 295780 

7(1,11= Pl( 1 1 * TAN ( ( ANG1 (II+ANGHI — 1) 1/2.) 

200 7.(1. 21= Pl( I)*TAN( I ANG1 ( IHANGICHI) 1/2.) 

Z (NSURF 1, 1 1 = PH NSURF 11*TAN( ( ANG1 (NSURFIHANGHNM) 1/2. I 
Z (NSURF 1, 21= 0. 

CALL TINTGR ( XI ,Z. NSURF 1 , SUM) 

F Xl= SUM) NSURF 1) 

J= 1 

WRI TF( 6.1000) 

WRITE) 6. 1250 ) J.J 
WRITE) 6.12901 

WRITE) 6, l 300 ) (Xl( I) ,P1( I ) , VI ( I ) . Z ( I ,11 , Z( I ,21 .SUM ( I ) , I = l , NS J RF1 ) 
00 210 1= 1. NSURF l 
7 ( I. 11= -PH II 
210 7.41,21= -PH I) 

CALL TINTGR) XI .Z.NSURF1, SUM) 

F Y 1= SUM(NSURFl) 

WRITE) 6,1260 ) J.J 
WRITE) 6.1290) 

WRITE) 6.1300 ) < XI ( I ) ,P1( I I ,V14 i) ,Z(I ,1) ,Z( 1 ,2) . SUM ( I ) , I = l , NS JRF 1 ) 
NM= NSURF 2- l 

IF (ANG2) D.EO.O.) ANG2( l ) =ANG S2( 1) 
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IF (ANG2! 2).E0.0.) ANG2(2 )=ANGS2<2 ) 

ANG2( 1.1= ANG 2! 1 1/57. 295780 
ANG2( 21= ANG2( 21/57.295780 
7(1.11= 0 . 

Z ( i* 2) = — P2 ( 1)*TAN( (ANG2(l)+ANG2(2))/2.) 

DO 220 1=2, NM 

IF 1 ANG2! I + ll.EO.O. ) ANG2 (I*1I=ANGS2(I«-1) 

ANG2(I*1> = ANG2( 1*11 /57. 295780 
7(1,11= -P2( II* TAN! ( ANG2 ( I ) +ANG2 ( I-ll 1/2.1 
220 7 (1.21= — P2( I ) *TAN( ( ANG2 ( I ) +ANG2 (I + l) )/2.) 

7 (NSURF 2, l) = -P2(MSURF2)*TAN( ( ANG2 { NSURF2)* ANG2< NM) ) /2. 1 
/ (NSURF2. 21= 0. 

CALI TINTGR! X2.Z, NSURF2. SUM) 

FX2= SUM( NSURF2 ) 

J= 2 

WRITE! 6, 1270 ) J, J 
WRITE! 6, 1290 ) 

WRITE! 6,1 100 ) ( X2(I) ,P2( I ) , V2( I) ,7(1 ,1) ,Z( I ,2) .SUM! I ) , 1= 1 , NS JRF2 ) 
DO 230 1=1 , NSURF2 
7(1,1)= P2( I ) 

2 30 7 (1,2)= P 2 ( I ) 

CALL T1NTGR(X2«Z , NSURF2 , SUM ) 

FY2= SUM(NSURF2I 
WRITE! 6, 1280 ) J,J 
WRITE! 6,1 290 ) 

WRI TE! 6.1300) (X2(I),P2(I)«V2( I) , Z! I , 1 ) ,2(1,2) ,SUM( I ), 1 = 1, NSJRF2) 
IF (KSURF.E0.2) GO TO 280 
NM= NSURF 3-1 

IF < ANG3! II ,E0. 0.) ANG3! 1 ) =ANG S3 ( 1 ) 

IF ( ANG 3! 2 ) . FQ. 0. I ANG3! 2 ) =ANG S3 ( 2 ) 

A NG3 ( l ) = ANG 3! 1J/57. 295780 
A NO 3! 2 ) = ANG3! 21/57.295780 

7(1.11= 0. 

Z ( 1 , 2) = P3( 1 )* TAN( ! ANG 3! 1 ) ♦ ANG 3(2) 1/2.1 
DO 240 1=2. NM 

IF ( AMG3! I + D.EO.O. ) ANG 3! 1*1) = ANG S3 ( I + l ) 

ANG3(H-1) = A NG3(I+l)/57. 295780 
7(1,1)= P3( I ) * TAN ( ( ANG3 ( I I+ANG3! 1 — 1) 1/2.1 

240 7(1,2)= P3( I)*TAN( (ANG3(I )+ANG3( I+l) )/2.) 

Z ( NSURF3, l)= P3(NSURF3)*TAN( (ANG 3 ( NSURF3 )♦ ANG3( NM) 1/2.1 
Z ( N SURF 3, 2 ) = 0. 

CALL T INTGR(X3,Z,NSURF3,SUM) 

F X3= SUM( NSURF 3 ) 

J = 3 

WRITE! 6.1250) J.J 
W R I TE( 6, 1290 ) 

'RITE! 6.1 300) (X3(I ) , P 3 ( I ),V3( 1) ,Z(I ,1) , Z ! I • 2 1 , SUM (I ) . I =1 , NSURF3 ) 
DO 250 1= 1.NSURF3 
2(1.1)= — P3( I) 

250 Z ( 1 , 21 = -P3( II 

CALL T INTGR ( X3 ,Z, NSURF 3* SUM) 

F Y3= SUM! NSURF3) 

WRITE! 6,1260) J.J 
WR I TF ( 6, 1 290 ) 

WRITE! 6, 13001 ( X3< I ) , P3 ( I ) , V3 (I) , Z ( I , 1 ) , Z ( I .2 ) .SUM ( I ) , I = 1 , NS J RF3 ) 
NM= NSURF 4-1 

IF ! ANG4! 1 ).F0.0.) ANG4! 1 ) =AN3 S4( 1 ) 

IF ( ANG4! 2).E0.0.) ANG4( 2 ) =ANG S4 ( 2 ) 
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ANG4I1)= ANG4! I ) /57. 295780 
ANG4! 2 ) = ANG4!2)/57. 295780 
Z ( 1, 1) = 0. 

Z 1 1.21 = -P4! 1 ) * TAN ( ( ANG4! 1 I+-ANG4! 2) )/2. » 

00 260 1=2. NM 

IF ( ANG4! I+II.EO.O.) ANG4I I ♦ 11 =ANGS4 C 1*1 ) 

ANUIU1I- ANG4I I«-l)/57.295780 
Z 11. 1> = -P4I 1 1 * TAN( ! ANG4! D+ANG4! 1-1) )/2.) 

260 Z ( 1 . 2 ) = -P4! I)*TANI IANG4I I H-ANG4! 0111/2.) 

Z INSURF4, 11= -P 4! N SURF 4) *TAN!!ANG4! NSURF4) + ANG4 CNM))/2.) 

Z I NSURF4. 2 ) = 0. 

CALL TINTGRI X4.Z.N SURF4.SUM) 

FX4= SUM! N SURF 4) 

J= 4 

WRITE! 6.1270) J.J 
WRITE! 6.1290 ) 

WRITE! 6, 1300) I X4( I ) , P4I I ) . V4t I ) . Z 1 1 , 1 ) , Z( I , 2 ) .SUM! I ) * 1 = 1. NSURF4) 
DO 270 1= 1.NSURF4 

1 I I. L)= P 4 ( I ) 

270 7 ( I.2)= P4I 1 ) 

CALL TINTGRI X4.Z.NSURF4, SUM) 

F Y4= . SUM! NSURF4) 

WRITE! 6,1280) J.J 
WRITE! 6, 1290 ) 

WRITE! 6. 1300) 1 X4( I ) .P4<1 ) , V4 1 1 ) .Z(I ,1) ,Z! I ,2) ,SU Ml I ) , I = 1 , NSJ RF4 ) 

ROTATE X AND Y FORCES TO OBTAIN COMPONENTS OF 
MERIDIONAL AND TANGENTIAL FORCES 

2 80 CALPH= CO SI A LPHMX ) 

S ALP H= SIN! ALPHMX) 

FMl= EX 1*CAL PH-FY l* SALPH 
F Tl= FXl* SALPH+FY 1*CALPH 
FM2= E X2*C AL PH— FY2* SAL PH 
F T2= FX2*SALPH«-FY2*CALPH 
IF (KS0RF.E0.2) GO TO 290 
FM3= FX3*CALPH-FY3*SALPH 
F T3= EX 3* SAL PH*-FY3*CALPH 
F M4= FX4*C ALPH— FY4*SALPH 
F T4= FX4* SALPH + FY4*CALPH 

ROTATE MERIDIONAL AND TANGENTIAL FORCES TO OBTAIN COMPONENTS QF 
LIFT AND DRAG FORCES 

290 VMI = VRI*COS(BETAI ) 

V T I = VR I* SI N !BE TA 1 ) 

V M0= VRO*COS(BETAQ) 

V T0= VRO* SIN(HETAQ) 

B ETAM = ATAN! ! VT I4-VT0 )/! VMI + VMU) ) *57. 295780 
WRI TEI 6, l 3 10 I BETAM 
BETAM= BETAM/57. 295780 
C BETA= CO SI BETAM) 

S 8ETA= SIN! 8ETAM) 

FDl= FMl*CBETA+FTl*SBETA 
EL 1 = -FMl*SB£TA*FTl*CaETA 
F 02= FM2*CBETA+FT2*SBETA 
FL2= -FM2*S8£TA+FT2*CBETA 
IF IKSORF.EO.2) WR I TE !6, 1330) 
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IF (KSURF.E0.4I Wft I TE ( 6* 1340) 

WRITE! 6,1320) FXI,FY1,FM1.FT1,FD1 .FL1 
IF (KSURF.E0.2) Wft I TE(6, 1350) 

IF (KSURF. EO. 4) WR I TE 4 6, 1 3601 
WR1TE( 6,13201 FX2.FY2.FM2.FT2,F02,FL2 
IF 4KSURF.EO.2I GO TO 300 
F 03= FM3*CBETA*FT3*SBETA 
F L3= -FM3*SBETA+FT3*CBETA 
F 04= FM4*CBETA+FT4*SBETA 
FL4= -FM4*SBETA+FT4*CBETA 
WRITE! 6.1370) 

WRITE! 6.1320) FX3.F Y3 , FM3 ,FT3 , FD3 , F L3 
WRITE (6, 1380) 

WRITE! 6, 1320) FX4. F Y4.FM4 .FT4, FD4 ,FL4 

SUM FORCES ON BLADE SEGMENTS, RATIO FORCES FOR TANOEM BLADES, 
AND COMPUTE TOTAL FORCES ON BLADE SECTION 

TWO SURFACES 
300 F Xl2 = F XI *F X 2 
FY12= FY1+FY2 
F Ml 2= FM l «-FM2 
F T12= FT1+FT2 
F012= F01+FU2 
FL 1 2= FL l «-FL 2 

IF (KSURF.E0.2) WR I TE ( 6, 1390) 

IF (KSURF. EO. 41 WR I TE !6. 1400) 

WRITE! 6,1320 ) FXI2.FY12.FM12 ,F T12.FDI 2 .FLI2 
TF12- SORT ( FM12**2*FT12**2I 
A TF 12M= ATAN(FT12/FM12)*57.295780 
IF (ATF12M.LT.0.) ATF12M= 180.+ATF12M 
IF ( KSURF .E0.2 ) Will TE (6,1410) TF12 
IF ( KSURF. E0.4) WR I TE ( 6, 1420 ) TF12 
WRITE! 6, 1430) ATF12M 
IF (KSURF. E0.4I GO TO 310 
GO TO 10 
C FOUR SURFACES 
310 FX34= FX3+FX4 
F Y34= FY3 +FY4 
FM34= FM3+FM4 
F T 34= FT3+FT4 
F 034= F03+FQ4 
FL34= FL3+FL4 
WRITE! 6,1440) 

W R I TF! 6,1320 ) FX34.F Y34.FM34.E T34 , FD34 .FL34 
TF34= S0RT(FM34**2+FT34**2) 

ATF34M= A TAN (FT34/FM34) *57.29 5 780 

IF (ATF34M.LT.0.) A TF34M= 180.+ATF34M 

WRITE! 6,1450) TF34 

WR ITE! 6. 1430 ) ATF34M 

F XR AT= FX 34/FX12 

F YR A T= FY34/FY 1 2 

FMRAT= FM34/FM 12 

F TRAT= FT 34/FT12 

FDRAT= F034/F012 

FLRAT= FL 34/ FL 1 2 

WRITE! 6,1460 ) F Xft AT , F YRA T ,FMRA T ,F TRAT , FDRAT , FLRAT 
TFRAT= TF 34/ TF 12 
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WRITE* 6, 1470) TFRAT 
F X 14— FX12+FX34 
F Y 14= Fy 1 2+F Y34 
FM14= FM12+FM34 
F T 14= F T1 2*-F T34 
FD14= F012+FD34 
FL 14= FL12+FL34 
WRITE* 6,1480) 

WRITE* 6.1320 ) FX14.F Y14.FM14.F T14 »FD1 4 »F LI 4 
T F14= S0RT*FM14**2+FT14**2) 

ATF14M = A TAN IF T14/FM14)*57. 295780 
IF IATF14M.LT.0.) ATF 14M= 180.+ATF14M 
WRITE* 6,1490 ) TF14 
WRITE* 6,1430) A TF14M 
GO TO 10 

320 WRI TE* 6,1 500 ) 

STOP 

FORMAT STATEMENTS 

1000 FORMAT* 1HI///) 

1010 FORMAT* 1615) 

1020 FORMAT* 4F 10.5) 

1030 FORMAT* 8F10.5) 

1040 FORMAT* 80H 

1 ) 

1050 FORMAT ( //6X,3HGAM,9X,1HR,11X, 3HPTZ ,8X,3HTTZ/3X,F7.3»3K,F9.2,3X» 

IF 10.2«2X,F9.2// ) 

1060 F0RMAT(6X, 3HVRI ,9X,3HVRG,8X,586ETAI , 7X , 5HBETA0, 7 X, 6HAL PHMX / 3X , 

IF 10.5* 2X.F10.5, 1X,F8.3,4X,F8.3»5X,F8.3//) 

1070 FORMAT* 6X , 3HKP V , 8X, 5HKSURF , 7X. 6HNSURF 1 , 6X .6HNSURF2 , 6X, 6HNS J RF3 , 6X, 
16HNSURF4/6X. 12 , l OX . I 2 , 10X . I 3 * 9 X, I 3 1 9 X , 1 3 ,9 X , I 3// ) 

1080 FORMAT* /5X, l 5HSUCTI0N SURFACE) 

1090 FORMAT* /5X, 16HPRESSURE SURFACE) 

1100 FORMAT* /5X.29HSUCTIQN SURFACE - FRONT BLADE) 

ULO FORMAT* /5X.3CHPRESSURE SURFACE - FRONT BLADE) 

1120 FORMAT* /5X.28HSUCTI0N SURFACE - REAR BLADE) 

1130 FORMAT* /5X.25HPRESSURE SURFACE - REAR BLADE) 

1140 FORMAT* 8X.2HXI , ILX.2HY1 ,HX,2hlPl ,10X,4HANGl) 

1150 FORMAT* 8X. 2HXI , 11X, 2HY l , 11X ,28 VI , 10X.4HANG1) 

1160 FORMAT* 8X.2HX2 . 11 X. 2HY2 , 1 1 X ,21 P2 , l OX , 4 HANG2 ) 

1170 FORMAT* 8X.2HX2. 11X,2HY2,11X,2HV2,10X, 4HANG2 ) 

1180 FORMAT* 8X.2HX3, 11 X.2HY3.1 IX, 2HP3 .10X.4HANG3) 

1190 FORMAT* 8X.2HX3, UX. 2HY3 , 1 1 X,2K V3 , 10X , 4HANG3 ) 

1200 FORMAT* 8X.2HX4, 11X.2HY4, 1 1 X.28P4 , 10X, 4HANG4 ) 

1210 FORMAT* 8X.2HX4, 1 1 X, 2H Y4 , 1 1 X ,28 V4 , 10X.4HANG4) 

12 20 FORMAT* 3X,F10.6,3X,F10.6,3X,F10.4,3X,F10.4) 

1230 FORMAT ( // /5X »8H SURF ACE ,11,208 - SPLINE FIT OF X-Y ,42X, 6HS PL I NE, 

19 X, 5HINPUT/8X, 1HX , l 1 X , 1H Y , 10X , 5HSLQPE ,7X , 1 1HSEC. DERIV..7X, 

29 HCURY ATURE ,6X, 5HANGLE, 10X.5H ANGLE ) 

12 40 FORMAT* * 3X , F 10. 6, 2X ,F 1 0. 6 , 1 X,E 12 . 5 , l X ,F1 5 .5 ,4X , FI 2 .5 , 3X , F9 . 4, 6X , 

1 F 9.4 ) ) 

1250 FORMAT* ///5X.8HSURFACE .11,658 - INTEGRATION OF {AVERAGE PRESSURE* 
1 T AN* AVERAGE ANGLE)) TO GET FX.Ii) 

1260 FORMAT* ✓/ / 5X . 8HSURF ACE ,11,328 - INTEGRATION OF * — P ) TO GET FY,I1) 
1270 FORMAT* ✓// 5X ,8H SURF ACE .11 ,668 - INTEGRATION OF ( -AVERAGE PRESSJRE 
l* TAN* AVERAGE ANGLE)) TO GET FX.II) 

1280 FORMAT* ✓/ / 5X .8HSU8F ACE .11,318 - INTEGRATION OF (P) TU GET FY,U) 
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l? 90 FORMAT ( 8X. 1HX. 10X,8HPRESSUR£, 6X.8H VELOCI T Y , 8X , 10 H I NTEGR AND 1 , 7X , 

1 1 OH IN TF GRAM 0 2.11 X.3H SUM) 

1300 FORMAT! { 3X . F 10. 6, 4X ,F 1 1. 5 ,4X , F 10. 3 ,3 X .F 14. 5 , 3 X ,F14 .5 , SX , FI 1 . 5 ) > 
1310 FORMAT! 1H 1// /4X, 61HANGLE BETWEEN MERIDIONAL-TANGENTIAL AXES AND LI 
IF T- DR AG AXES =.F9. A///////) 

1320 FORMAT! 12 X. lhX . 14X . 1HY , 9X , 1 OHMERI Dt QNAL , 5X » 1 OHT ANGENT I AL , 8X , 4HDRA3 
1. 11X.4HLIFT/6! 10X.5HFORCE ) /1X,6(4X,F11.4)//) 

1330 FORMAT! 4X.22FiSUCTIUN SURFACE FORCES) 

1340 FORMAT! 4X.36HSUCTI0N SURFACE FORCES - FRONT BLADE) 

1350 FORMAT! 4X . 23HPRE SSURE SURFACE FORCES) 

1360 FORMAT ! 4X • 3 7 HP RE SSURE SURFACE FORCES - FRONT BLADE) 

1370 FORMAT! 4X *35 HSUCTION SURFACE FURCES - REAR SLADE) 

1380 F0RMAT14X.36HPRESSURE SURFACE FORCES - REAR BLADE) 

1390 FORMAT! ////74X.50HSUM OF SUCTION SURFACE AND PRESSURE SURFACE FQRC 
IF S ) 

1400 FORMAT! /////4X.64HSUM OF SUCTION SURFACE AND PRESSURE SJRFACE FORC 
IE S - FRONT BLADE) 

1410 FORMAT! 9X.22FT0TAL FORCE ON BLADE =,F11.4) 

1420 FORMAT! 9X.28HT0TAL FORCE ON FRONT BLADE =,Fll.4) 

1430 FORMAT! 9X.39HANGLE WITH RESPECT TO MERIDIONAL AXIS =,Fll. W//) 
1440 FORMAT! 4X.63HSUM OF SUCTION SURFACE AND PRESSURE SURFACE FORCES - 
1R EAR BLADE) 

1450 FORMAT! 9X.27HT0TAL FORCE ON REAR BLADE =,F11.4) 

1460 FORMAT! 1H1///4X.48HRATIQ OF REAR BLADE FORCES TO FRONT BLADE FORCE 
1S/12X. 1HX » 14 X» IHY, 9X , 10HMER10 I ONAL . 5 X » 1QHT ANGE NTIAL*8X,4NDRAG» 1 IX, 
24HL IFT/61 9X , 6HF0RCE S) / 1 X.6 i 4X.F1 1 . 4) //) 

1470 FORMAT i 9X » 66FR A T 1 0 OF TOTAL FORCE ON REAR BLADE TO TOTAL FORCE ON 
1 FRONT BLADE =.F10.4 ///////) 

1480 FORMAT! 4X . 38FSUM OF ALL FRONT AND REAR BLADE FORCES) 

1490 FORMAT! 9X.36HT0TAL FORCE ON ENTIRE TANDEM BLADE =,Fll.4> 

1500 F OR MAT ! // / 1 X .29HERR0R IN INPUT - KPV OR KSURF ) 

END 


SlBFTC SPUN 

SUBROUTINE SPL INFIX. Y,N,DY0X,D2 YDX2 .CUR V * ANG) 

SPLINE FITS A SPLINE CURVE TO X AND Y, AND CALCULATES FIRST AND 
SECOND DERIVATIVES, CURVATURES, AND ANGLES AT THE SPLINE POINTS 

D IMENS 1 ON X!N> .YIN) .OYDX(N) ,D2YOX2(N> .CURV(N) , ANG! N) 

DIMENSION GI 100 ) »H ! 100) 

G! 1 )= -2.0 
H! 11= O. 

N l— N-l 

IF INI.LT.2) GO TO 20 
00 10 1=2 »N 1 
A = ( XI I )-X( I — 1 I ) / 6. 

B= < XI I-*l )-X!I ) )/6. 

C = 2.*! A+B I- A*G 1 1-1 ) 

0= 1 Y1 1*1 )-YU ) )/!Xl 1+D-X1I ) 1-IY1 I)-Y1 1-1) )/( X( I )-X( i-1) ) 

G 1 1 1 = B /C 

10 H ! I ) = ! 0- A*H ( I — 1 ) ) /C 
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20 0 2YQX21 N ) -= H (N 1 ) / 1 • 5+G (Nil) 

00 30 1=2, N 
K = N+l-I 

30 D2YDX?<K) = H (K ) — G < K)*02YQX24K* II 

0 YOX ( 1»= i.X ( II — XC 2» I /6.* ( 2^*02 Y0X2 11) *02 YDX2 ( 2 ) ) ♦( Y ( 2 ) -Y < 1) ) / { X< 2 ) 
1- X ( 1 ) ) 

00 40 1=2. N 

40 0 YOX ( I ) = ( X ( 1)- XI 1-1) ) /6.*<2.*02YDX2( I ) +02 YOX 2 (1-1 ))+(Y(I)-T( 1 — 1)1 
1/ < X ( 1 1 — XI 1-1 )) 

00 50 1=1, N 

CURV(f) = 02Y0X2UI /< L.+OYOXm **2 I **1 • 5 
50 A NG( I ) = A rANlOYDX(I) >*57.295780 
RETURN 
END 


$ I BE TC T INTO 

SUBROUTINE T IN TGR ( X , Y ,N , SUM) 

TINTGR INTEGRATES THE SET Uf DATA POINTS X-Y USING THE 
TRAP AZO IOAL RULE 

0 IMENSION X( 100) ,Y( 100,2) .SUM! 100) 

S IJM ( 1)= 0. 

00 10 1=2, N 

10 S UM ( I ) = SUM! I-D + IYI I » 1 )♦ Y( I— l ,2) ) *(XII)— X( l-l ) )/2. 
RETURN 
END 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, July 7, 1970, 

126-15. 
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APPENDIX -SYMBOLS 


F force on blade section or blg.de 

surface, per unit of blade 
height, lbf/ft; N/m 

M distance in meridional direc- 

tion, ft; m 

P pressure, lbf/ft^; N/m^ 

R distance in radial direction, 

ft; m 

T temperature, °R; K 

V relative velocity (figs. 3 

and 6), ft/sec; m/sec 

mean relative velocity across 

ixiG3tn 

blade section (fig. 6), 
ft/sec; m/sec 

X distance in X (input) direc- 

tion, ft; m 

Y distance in Y (input) direc- 

tion, ft; m 

Z distance in axial direction, 

ft; m 

a blade surface angle with re- 

spect to the X axis (fig. 5), 
deg 

«Mx angle between meridional axis 

and X axis of input (fig. 2), 
deg 

/3 relative flow angle with re- 

spect to meridional or axial 
direction (figs. 3 and 6), deg 

^mean mean flow angle across blade 
section (fig. 6), deg 


y ratio of specific heats 

9 distance in tangential direc- 

tion, rad 

o o 

p density, slug/ft ; kg/m 

Subscripts: 

D in the drag force direction 

(parallel to direction of mean 
velocity) 

i at a local point 

in at the blade inlet 

L in the lift force direction (per- 

pendicular to direction of 
mean velocity) 

M in the meridional direction 

out at the blade outlet 

p pressure surface value 

s suction surface value 

tot total force 

X in the X (input) direction 

Y in the Y (input) direction 

Z in the axial direction 

9 in the tangential direction 

1, 2, 3, 4 identification of blade sur- 
faces 1, 2, 3, or 4 

Superscript: 

' total or stagnation value 
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